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ABSTRACT 

In this work we present the first results of our imaging campaign at Keck Observatory to identify 
the host galaxies of "dark" gamma-ray bursts (GRBs), events with no detected optical afterglow 
or with detected optical flux significantly fainter than expected from the observed X-ray afterglow. 
We find that out of a uniform sample of 29 Swift bursts rapidly observed by the Palomar 60-inch 
telescope through March 2008 (14 of which we classify as dark), all events have either a detected 
optical afterglow, a probable optical host-galaxy detection, or both. Our results constrain the fraction 
of Swift GRBs coming from very high redshift (z > 7), such as the recent GRB 090423, to between 0.2- 
7 percent at 80% confidence. In contrast, a significant fraction of the sample requires large extinction 
columns (host-frame Ay ^ 1 mag, with several events showing Ay > 2 — 6 mag), identifying dust 
extinction as the dominant cause of the dark GRB phenomenon. We infer that a significant fraction of 
GRBs (and, by association, of high-mass star formation) occurs in highly obscured regions. However, 
the host galaxies of dark GRBs seem to have normal optical colors, suggesting that the source of 
obscuring dust is local to the vicinity of the GRB progenitor or highly unevenly distributed within 
the host galaxy. 

Subject headings: gamma-rays: bursts — dust, extinction — galaxies: photometry — galaxies: high- 
redshift 



1. INTRODUCTION 

The prevalence of "dark" gamma-ray bursts (GRBs) 
remains one of the most persistent mysteries of the 
field, twelve years a f ter th e discovery of GRB afterglows 
(jvan Paradijs et all 119971 : ICosta et al.l [1993). While 
we now know that GRBs are frequently accompanied 
by extremely lu minous afterg lows (some times spectac- 

!-g-, lAk a 

iRacusin et al.l[2008f l an optical detection is reported in 
only about half of cases since the launch of Swift 
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ularly so: e.g lAk erlofejalJ [l999;' Bl oom et all l2009l : 
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In contrast, an X-ra y detection is nearly always re- 
ported for Swi/t bursts (|Gehrelsl | 2008l ). Partly this is due 
to observational constraints: the limitations of ground- 
based observing prevent a significant fraction of GRBs 
from being observed with terrestrial optical telescopes 
at all. F urthermore, the Ult ra- Violet /Optical Telescope 
(UVOT: 'Romi ng et al.l[2005l ) on-board Swift has a typi- 
cal limiting magnitude that i s shallower than the equiva- 
lent X-Ray Telescope (XRT; iBurrows et "alll2005[) X-ray 
fiux limit for a typical broadband afterglow spectrum, 
in particular when filters are applied. Galactic extinc- 
tion, stellar crowding, and proximity to the Sun or Moon, 
which do not significantly affect the X-ray band, also of- 
ten complicate optical follow-up. Estimates for the in- 
trinsic frequency of optically dim GRBs vary and likely 
depend on the sensitivity of the detect ing satellite, but 
for Swift events lAkerlof fc SwanI (|2007| ) have estimated 
that approximately 30% of GRBs have an optical mag- 
nitude > 22 at only 1000 s after the trigger, and 15-20% 
have an optical magnitude > 24 at this time. Detecting 
an optical afterglow from such an event requires a rapid 
response by a large-aperture telescope and is rare. 

It is noteworthy that most of the conclusions about 
GRBs to date are based on a limited subsample of 
well-studied events that tends to exclude this large 
population of faint afterglows. For example, evidence 
of a GRB-SN connection can only be established for 
known low-redshift events targeted for inte nsive photo- 
metric and spectroscopic follow-up (but c.f. iLevan et al.l 
i200 5). Likewise, conclu sions based on the nature o f GRB 
host galaxies ( Bloom c t al.l 120021 : iFruchter et all l200d 
IWainwright et ah. 2007.) require accurate (generally sub- 
1") positions. Only a handful of pre-Swift events without 
optical counterparts had sufficiently precise positions for 
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later follow-up work of this nature. Therefore the specific 
scrutiny of optically dark events is vital to understanding 
the entire GRB demography. 

Key in the study of dark bursts has been the 
progression from a wholly observational definition of 
darkness to the physic ally motivated f3ox criterion of 
iJakobsson et al.l (|2004af) . who define a dark burst on the 
basis of the flux ratio between X-ray and optical band- 
passes in the afterglow at 11 hr after the burst. Here the 
parameter /3ox is the observed spectral index (defined 
using the convention oc between the X-ray and 
optical bands, after correcting for Galactic extinction: 
Pox = log(Fx/Fopt)/log(Ax/Aopt)- Jakobsson defines a 
dark burst as one with /3ox < 0.5, motivated by the pre- 
diction from the synchrotron model in which, once the 
afterglow begins to fade, the intrinsic spectrum is given 
by F oc (for < i^c) or F (x v~p/'^ (for v > Vc), 

implying /3ox > 0.5 if p > 2. 

The availability and uniformity of X-ray follow-up in 
the Swift era makes this definition of darkness particu- 
larly appropriate for a survey of Swift bursts. Even so, 
a purely optically defined criterion is still relevant: opti- 
cal detection versus non-detection (rather than the flux 
ratio) is an essential factor determining the nature of fur- 
ther follow-up of the event: sensitive searches for host- 
galaxy dust, spectroscopic redshifts and measurements 
of the host ISM properties, and (to a lesser extent) accu- 
rate host identification require bright optical afterglows, 
making this likely the dominant selection bias affecting 
our current understanding of GRB afterglows and their 
origins. 

The implications of dark bursts are potentially far- 
reaching, and the importance of folding them into our 
understanding of the GRB population as a whole is great, 
as — depending on the cause(s) of their optical faintness — 
there are reasons to suspect that their nature or envi- 
ronments may differ from those of the optically brighter 
GRBs which underpin our underst anding of the field . 
Some of the possibilities include fe.g.. lFvnbo et al1l200lD : 

1. Extinction. Dust in the GRB host galaxy (or 
elsewhere along the line of sight) can strongly 
obscure the rest-frame optical and ultraviolet 
light, dimming and redde i iing the afterglow 
fe.g.. iDiorgovski et alJ 120011: iLazzati et all l2002t 
iReichart fc Pried l2002f ). While previous (largely 
optically selected) samples have shown little ev- 
idcnc e for widespread dust along GRB sightline s 
(e.g., IKann et aLl[2006l . 120071: ISchadv et alJ[200l . 
recent cases such as GRB 080607 {Ay = 3.2 
mag^^, iProchaska et al.l [20091 ) have demonstrated 
that very large dust columns can and do occur. 
A bias against dusty galaxies in the current sample 
could easily mislead us in conclusions about, for ex- 
ampl e, mean GRB host me t allicities and luminosi- 
ties (jFruchter et all 120061: iWolf fc Podsiadlowskil 
[2007h . 

2. High redshift. GRBs have now been observed 

In addition to assuming p > 2, this definition is meaningful 
only if the synchrotron model is assumed to be a complete descrip- 
tion of the afterglow SED at these wavelengths. We will make these 
assumptions throughout the paper. 

Throughout this paper, Ay refers to extinction in the host 
galaxy rest-frame V-band. 



out to z = 8.3 ()Tanvir et al.ll2009l : ISalvaterra et al.l 
[2009). At z > 6, photons which would be red- 
shifted into the optical bandpass are absorbed by 
neutral hydrogen in the host galaxy and IGM, sup- 
pressing th e observed opti cal flux almost entirely 
()Gunn fc Peterson 19 65: Fan et al.l[2006,) . The red- 
shift distribution of GRBs beyond z ~ 6 (and 
its implications on the star-formation history of 
the universe) cannot be observationally constrained 
without incorporating the dark burst population. 

3. Low luminosity. It is well-esta blished (e.g., 
iGehrels et all 120081 : lNvsewander~eral, ,2009) that 
GRB fluence and afterglow flux are positively cor- 
related (that is, underluminous bursts tend to also 
have underluminous afterglows) . Due to a wide dis- 
tribution both in the depth of optical follow-up as 
well as in the gamma-ray fluence of observed GRBs, 
many nondetections could simply be attributed to 
follow-up that was not deep enough to constrain 
the predicted optical afterglow for a relatively faint 
GRB, without need to invoke absorption effects. 

4. Low-density medium. However, it is physically 
possible to have a energetic event without a lumi- 
nous afterglow. The afterglow phenomenon, which 
is thought to originate from shocks in the surround- 
ing medium (jPaczvnski fc Rhoadslll993l ). critically 
depends on the presence of circumstellar gas at suf- 
ficient density to excite bright synchrotron radia- 
tion. GRBs exploding in galaxy halos or the inter- 
galactic medium are predicted to have afterglows 
orders of magnitude fainter than those occurrin g 
in galactic disks fe.g. jKumar fc Panaitescull2000D . 

To some extent, these various possibilities can be dis- 
entangled via broadband observations of the afterglows 
of the events alone. For example, a low-density medium 
will result in a dim afterglow at all wavelengths, extinc- 
tion will suppress both the optical and the near-IR flux 
as well as soft X-rays (to different and characteristic ex- 
tents), and a high redshift will suppress only the optical 
flux. As a result, we will give attention in the subse- 
quent discussion to the nature of the afterglows at all 
wavelengths. However, extensive broadband follow-up is 
not always available (and the decision to trigger multi- 
wavelength observations carries its own selection biases), 
and in some cases the two possibilities are difficult to 
disentangle. 

The remaining degeneracies can largely be broken via 
deep imaging of the host galaxy of a GRB. In particular, 
high-redshift bursts should not have optically observable 
host galaxies, and the detection of a host can rule out 
the high-redshift hypothesis for that event. Secondar- 
ily, study of the host galaxies themselves can determine 
whether our existing sample of pie-Swift host galaxies is 
in fact typical, or if we are missing (for example) a large 
population of red, dusty ULIRGs. 

2. THE PALOMAR 60-INCH SAMPLE 

The Palomar 60-inch telescope (P60; ICenko et al.l 
2006a) is a robotic facility designed for moderately fast 
{t < 3 minutes) and sustained (i? < 23 mag) observations 
of GRB afterglows and other transient events. Fully op- 
erational since 2004 September, the P60 now routinely 
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interrupts regular queue- scheduled observations in re- 
sponse to electronic notification of transient events. The 
standard P60 response to Swift GRB alerts results in a 
sequence of multi-color (gRci'z') observations for ap- 
proximately the first hour after the trigger. Subsequent 
observations are then triggered manually based on the 
properties of the afterglow in observations to that point. 

The firs t catalog of P60 G RB observations was pre- 
sented bv iCenko et"aLl (|2009D . The P60 follow-up pro- 
gram is fully robotic and the GRBs presented in that 
sample were selected entirely based on whether an event 
was rapidly foUowed-up. P60 automatically follows up 
all Swift GRB triggers that are observable, and therefore 
this catalog constitutes an effectively uniform sample of 
Swift events to date, and should not be affected by any 
afterglow-related biases. Other advantages of this popu- 
lation include a high afterglow detection efficiency (75%, 
thanks to the relatively large aperture of the telescope 
and red filter sequence) and a large fraction with spec- 
troscopic redshifts (60%). 

In total, the P60 sample contains 29 events (Table [T]). 
Of these, 7 were undetected with the P60 (to a typical 
limiting magnitude of i? > 20 — 23, depending on con- 
ditions) at 1000 seconds. No event that was undetected 
at 1000 s was detected at earlier times. This is approx- 
imately consistent with the results of previous studies 
which have attempted to correct for the shallow follow- 
up of most Swift GRBs in determining t he true afterglow 
bright ness function: in particular that of lAkerlof fc Swanl 
()2007D . which estimates (Figure 6 of that work) that 30% 
of afterglows are fainter than 22nd magnitude at 1000 s. 
These events are "dark" by the simple nondetection cri- 
terion, although the rapid response, large aperture, and 
nearly uniform depth of P60 makes a nondetection signif- 
icantly more meaningful than is typical for Swift bursts 
(many of which have no optical follow-up at all, or follow- 
up only from the UVOT and small-aperture ground- 
based telescopes.) Four of the seven events have optical 
or infrared afterglows detected by other telescopes (typ- 
ically with larger apertures and/or a redder wavelength 
response.) 

We include a handful of additional eve nts as "dark" via 
apphca tion of the /3ox < 0.5 criterion of lJakobsson et all 
((2Q043), though we apply it at 1000 s instead of 11 hr, 
given that late-time imaging is not always available and 
that our non-detection cutoff is also at 1000 s. There 
are 12 such events that satisfy this criterion: 5 of which 
are also P60 nondetections and 7 events which are de- 
tected by P60, but at a flux level that is less than a 
simple (3 = 0.5 extrapola tion of the 2 keV X -ray flux as 
determined by Table 3 in ICenko etall (|2009D . There- 
fore our full "dark" sample defined by the union of both 

This involves some risks: there are occasional cases in which 
X-ray rebrightenings or strong spectral evolution is observed af- 
ter 1000 s, indicating the contribution of additional prompt-like 
emission (X-ray fl ares) which have m uch harder spectra than a 
typical afterglow I IButler fc Kocevskil ]2007b) and could generate 
"pseudo-dark" events at early times which would look normal in 
later observations. We will discuss the possibility of this contribu- 
tion in the next section in the few cases where there appears to be 
evidence of extended activity at this time, but conclude that it is 
not a significant contaminant of our dark burst sample. 

T wo events are listed with /3ox < 0.5 at 1000s in'Ce nko et al.l 
II2009I') which we do not include in our sample: GRB 050820A and 
GRB 071003. In both cases, the Swift XRT was not observing the 
source at 1000s and the actual spectral index at that time is un- 



criteria consists of 14 events in all, approximately half of 
the P60 sample. All 14 fields were imaged to deep limits 
at Keck Observatory, as discussed in the next section. 

3. OBSERVATIONS 

3.1. The Keck Imaging Campaign 

Since 2005, we have been acquiring d eep optical imag,- 
ing of GRB fields using LRIS on Keck I (lOke et allUQoH) 
as part of our ongoing Berkeley Keck GRB Host Project. 
Primary goals of the survey include elucidating the ori- 
gins of dark GR Bs, studying the hosts of X-ray fiashes 
(XRFs; iHeise e t al. 2001) and short gamma-ray bursts 
(SHBs), constraining late-time supernova emission from 
bursts, and i dentifying the hos ts of GRBs showing strong 
DLAs (e. g., iChen et al.l 120091) or intervening Mg H ab- 
sorption (i Pollack et al.l [20091 ) . As part of this project, 
we conducted imaging of all 14 P60 "dark" bursts above. 
These observations were supplemented in a few cases by 
additional imaging taken by the Caltech GRB group, also 
with LRIS. 

Our typical imaging mode was with the R and g fil- 
ters simultaneously with the D560 dichroic, but other 
setups were also frequently used, with varying expo- 
sure times and total integrations depending on the 
field. Observations were conducted between 2005 and 
2009 across 15 different observing runs. Some of these 
nights were photometric, and the photometry was cal- 
ibrated using Landolt standard fields (SA 92, SA 101, 
PG 2213, and Markarian A: iLandoltl Il992f ). Non- 
photometric nights were calibrated using Sloan Digital 
Sky Survey (SDSS) phot ometry of stars in the GRB field 
(|Adelman-McCarthv et a l. 2008), or when unavailable, 
using cali brations from th e P60 matched to USNO stan- 
dards (Monet et al.ll200ai . One field (GRB 060210) was 
calibrated to USNO directly {B2, i?2, and I magnitudes), 
using an average color oi 13 — g ^ 0.4 ± 0.3 (jJordi et al.l 
120061 ) to convert B to g. A list of all observations and 
exposure times is presented in Table [2] 

Images were reduced using standard techniques using a 
custom Pyraf script written in Python. Astrometry was 
conducted relative to USNO-Bl.O astrometric standards. 
In cases where we detected the optical afterglow with 
P60 (and, in on e case, with the robotic infrared telescope 
PAIRITEL: Bloom et al.ll2006f) these early-time images 
were registered and aligned with the Keck data to deter- 
mine the most accurate possible afterglow position rela- 
tive to the host-galaxy candidates. 

3.2. Host Identification 

Until recently, the same biases that made pre- 
Swift host searches difficult without optical positions 
have applied to Swift as well: early XRT positions were 
accurate to only 4-6" , an error region sufficiently large as 
to normally contain numerous faint galaxi es. However , 
by using optical sources to register the field ()Butlerll2007l : 

known; the estimate in lCenko et al.l l|2009l ) was based on an extrap- 
olation from other epochs. This is difficult, since GRB 050820A 
shows extensive early-time X-ray flaring while GRB 071003 experi- 
ences a dramatic rebrightening at around 1 day when XRT obser- 
vations begin. Late-time observations in both cases (Ccnko ct al] 
12006b : iPerlev et al.l l2009al) show that the spectral index is quite 
normal at late times, strongly indicating that neither event is a 
genuine dark burst by either of our criteria (these are, in fact, 
among the two brightest bursts of the Swift era. ) 
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TABLE 1 
P60 GRBs 



GRB Fluence^' X-Ray Flux'' R''''^ NIR^'d,c ^^^f Reason for dark 

(10~^ erg cm~^) (^Jy) (ma^g) (m'^'g) classification 



Dark GRBs 














050412 


6.18 


0.27 


> 21.4 




< 0.49 


P60 nondetection 


050416A 


3.67 


1.95 


20.31 




0.35 


low /3ox 


050607 


5.92 


0.45 


~ 22.1' 




~ 0.33 


P60 nondetection 


050713A 


51.1 


14.51 


18.45 




0.31 


low Sr^v 


050915A 


8.5 


0.72 


> 20.7 


~ 18 


< 0.44 


P60 nondetection 


060210 


76.6 


12.23 


18.2 




0.37 


low rJr-iv 


nfin^i OR 

U\JU0 L\JD 




1 ^ riQ 


~ 20 4^ 






low /3ox 




0.72 


0.17 


> 19.9 




< 0.76 


JT litJlivlt: LCljLrltJli 




S fiQ 
o.uy 


n Q9 


\ 99 
^ zz.u 


K ~ 18^ 


^ n 94 


P60 iiondctGctrion 


Ofil 999 A 


7Q Q 




> 22 1 


K IS 


^ n 1 Q"^ 
^ — u.iy 


P60 iiondctGctrion 


(Tync^o 1 
U/Uozl 


80.1 


4.40 


> 22.9^ 


K > 18.7 


< —0.10 


P60 nondetection 


noriQ 1 n A 

UoUolyA 


48 


1.19 


20.43 




0.41 


low /3ox 


noriQ 1 c\i^ 

UoUolyO 


36 


11.68 


18.32 




0.36 


low ^ox 


UsUoiU 


2.7 


1.37 


> 21.0 


z = 20.0 


< 0.31 


low ^ox 


Other GRBs 














050820A 


34.4 


~ 150S 


15.21 




~0.4 




050908 


0.51 


0.12 


19.17 




0.91 




060110 


15.7 


7.42 


15.46 




0.80 




060502A 


23.1 


1.22 


19.50 




0.53 




060906 


22.1 


0.20 


18.84 




0.88 




060908 


28.0 


0.92 


17.59 




0.82 




070208 


47.7 


0.88 


19.74 




0.54 




070419A 


5.58 


0.17 


19.02 




0.87 




071003 


83 


_h 


17.06 




_h 




071010A 


2.0 


2.11 


16.18 




0.89 




071011 


0.22 


8.06 


16.42 




0.66 




071020 


23 


6.91 


17.66 




0.52 




071122 


5.8 


0.34 


20.02 




0.64 




080310 


23 


2.19 


16.88 




0.79 




080319B 


810 


265.8 


13.69 




0.52 





low /3ox 



low /3ox 
low /3ox 
low /3ox 



15 - 150 keV; taken from the BAT GRB table. 

Absorbed flux at 2 keV; calculated using the Swift XRT Repository dEvans et al.ll2007l1 . 
" Calculated at 1000 s. 

Vega mag; corrected for Galactic extinction. 
" Specified only in the case of i?-band nondetections. 

^ Between ij-band and 2 keV. From Cenko et al. (2009), modified include deeper non-P60 upper limits (where available) and revised 
XRT light curves. 

s The XRT was not observing this burst at 1000 s, and earlier observations were dominated by rapid flaring (see footnote in text). 
^ The XRT did not slew to this burst until 22000 sec after the BAT trigger (see footnote in text). 
' IRhoads et~all ll2005l ') 

J Based on extrapolation from later times: the burst was not detected in i?-band at 1000s. 
"^ ITanvir et"all H^OB ) 

' Interpolated b e tween P 60 measurements and lRau et al.l 1120071) . 
"' ICenko fcTtoJ Il2006bl) 



|Go ad _ et al.| [2007HEvans et al.ll2009D . the Swift XRT now 
routinely produces afterglow positions to better than 2" 
(90% confidence). Furthermore, thanks to the prolifer- 
ation of small- to medium-sized telescopes and the im- 
proving ability of larger apertures to respond relatively 
quickly, all but three of the P60-followed bursts in our 
sample are detected in the optical or IR. In all cases 
where a host candidate is identified in or near the error 
circle, we follow the prescription in iBloom et al.l ()2002f ) 
to estimate Pchanco- Formally, this parameter is an esti- 
mate of the probability that one or more galaxies with 
an observed magnitude brighter than m will be centered 
within a randomly chosen region on the sky with solid 
angle tt6'^. This probability is given by: 

f ch = 1 - exp(7r6'^(7<,„) 

Where a<m is the average sky surface density of galax- 
ies with apparent magnitude brighter than to, taken in 



this case from iHogg et all ()1997[ ). The values for to 
and 9 for ea c h bur st-host association are chosen as in 
iBloom et al.l ()2002D . with two exceptions. Because we 
do not have access to space-based imaging and the size 
of a typical host galaxy is significantly smaller than the 
seeing disk, we conservatively use the visible extent of the 
optical disk in the ground-based imaging rather than the 
half-hght radius. We also use the 90% confidence radius, 
rather than 3cr, which is slightly less conservative. We 
treat this value as an estimate of the probability that, 
for a given burst, the association with the nearest host 
galaxy is incorrect. 

Some additional caution is warranted before interpret- 
ing Pchanco this way. In particular, this probability ap- 
plies to a single event treated in isolation only: it is not 
necessarily appropriate for events chosen from a larger 
sample which includes both detections and nondetections 
(a shallow survey of a very large number of well-localized 
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TABLE 2 
Keck Imaging Observations 



OHri t leld 


(Jbs. Date 
(U 1) 


Filter 


T 

Int. 


Seeing 
( ) 


Cal. Sys. 


Cal. Unc. 

(mag) 


5(T Limit 

(mag) 


Eb-v 
(mag) 


050412 


2007-12-13 


g 


690 


1.4 


SDSS'^ 


0.03 


25.8 


0.02 






R 


600 


1.2 




0.17 


24.5 




050416A 


2005-06-05 


g 


960 


0.9 


SDSS 


0.03 


26.2 


0.03 






R 


960 


0.9 




0.02 


25.4 




050607 


2007-10-09 





960 


1.0 


Landolf* 


0.3s 


24.4 


0.156 






R 


870 


1.0 




0.3 


23.7 




050713A 


2008-08-02 


(1 

ij 


990 


0.8 


P60/USNO'' 


0.25 


25.7 


0.414 






R 


870 


0.7 




0.27 


24.7 




050915A 


2005-12-04 


V 


2280 


0.7 


Landolt 


0.05 


25.8 


0.026 






I 


1539 


0.8 




0.02 


24.9 






2005-10-31 





1680 


1.0 


P60/USNO 


0.25 


25.5 








R 


1500 


1.0 




0.35 


24.5 




06021 


2007-08-1 3 


R 


540 


0.7 


USNO' 


0.35 


23.6 


0.093 




9nnQ 09 1 Q 


9 


1680 


0.8 




35 


24.4 








I 


1530 


1.0 




0.14 


23.5 




UUUOJ-U-D 


ZUUU-UO-O-i- 


9 


oo^u 


1 4 


ijandolt 


u.uz 


zo.o 


u.uoy 






R 


3660 


1.4 




0.02 


25.5 






2008-09-1 9 


9 


1080 


1.0 


SDSS 


0.04 


26.3 


0.024 






R 


1260 


1.0 




0.10 


24.8 




0fi0Q2'?A 


2007-04-1 fi 


Y 


1560 


1.4 


SDSS 


0.04 


25.2 


0.060 






J 


1590 


1.2 




06 


23.8 






9007 08 1 9 




1500 


0.8 




0.07 


26.4 










1500 


0.6 




0.09 


23.6 






9007 07 1 8 


y 


710 


0.8 




05 


24.7 


0.099 






I 


600 


0.7 




0.05 


23.7 






9007 08 1 9 




1500 


0.7 


PfiO /TTSNO 


0.12 


25.9 










1500 


0.6 




27 


23.6 






2009-05-31 


H 


900 


0.5 


2MASS 


0.06 


21.6 








K 


1800 


0.5 




0.09 


21.7 




070521 


2007-05-21 


V 


1500 


0.7 


SDSS 


0.05 


24.8 


0.027 






I 


1500 


0.8 




0.03 


24.3 






2009-06-25 


V 


1440 


0.7 


SDSS 


0.05 


26.2 








RG850 


1260 


0.8 




0.15 


24.6 




080319A 


2009-02-19 


9 


1070 


0.6 


SDSS 


0.07 


26.4 


0.015 






R 


960 


0.7 




0.04 


25.0 




080319C 


2009-02-19 


9 


1530 


0.9 


SDSS 


0.05 


25.6 


0.026 






R 


1380 


0.7 




0.13 


24.5 




080320 


2009-02-19 


9 


990 


1.0 


SDSS 


0.18 


25.8 


0.014 






I 


810 


1.3 




0.09 


24.1 





Total integration time. 

As measured over a 1" aperture and averaged over the field; not corrected for extinction. BVRI 
magnitudes are in the Vega system. The JIG850 filter i s calib rated to the SDSS 2-band. 
" Sloan Digital Sky Survey: lAdelman-McCarthv eFaLl ll200g) 

LandoTl 1199^) 

° P60 caUbration, based on USNO Bl.O catalog HMonet et al.ll200^ ') 
' Direct calibration to USNO Bl.O catalog. 

s The two standard star observations during the 2007-10-09 run are not consistent with each other, indi- 
cating that this night may not have been photometric. 



objects would find many individual low-Pchancc galaxies 
even if the positions were chosen completely randomly). 
Fortunately, in our case we identify good host galaxy 
candidates for most of our objects: 11 out of 14 fields 
contain at least one object with Pchanco < 0.1 consistent 
with the error circle. Nevertheless, given the number 
of fields observed, we must recognize that the chance of 
a misidentification being present somewhere in the full 
sample is not insignificant. A basic Monte Carlo analysis 
(including the nondetections) suggests that the probabil- 
ity of at least one chance coincidence being present in our 
host sample is an appreciable 48%, and the probability 
of two or more is about 15%. 

The Pchance Calculation also assumes that lines of sight 
toward GRBs, and in particular toward dark GRBs, 
are randomly sampled among all sightlines in the uni- 
verse. One possible interpretation of the overabun- 



dance of Mg II absorbers in GRB spectra relative to 
QSOs (jProchter et al.l I2006D is that this assumption 
is incorrect and observed GRBs preferentially cluster 
along lines of sight near low-2 galaxies, perhaps due 
to gravitational lensin g. This interpretation is gener- 
ally disfavored ( ,Prochter et al.l 1200(1 ) , and for the few 
cases of gal axy-associated Mg II systems in GRB spec- 
tra to date (iMase tti et"alll2003t iJakobsson et aI1l2004bl : 
iPollack et a l. 2009) there has been no clear demonstra- 
tion that the number and offset distribution of these 
galaxies implies a significant excess of what is expected 
from chance. Another possibility which could affect our 
results is if dark GRBs are due to extinction in unre- 
lated field galaxies along the light of sight (rather than 
in the host galaxy) and therefore more likely to fall close 
to a line-of-sight galaxy: such an effect was studied as a 
possible interpretation of the GRB-QSO discrepancy in 
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terms of a selection bias (|Sudilovskv et al.ll2007D . Were 
this the case, dark GRB sight lines would be biased to- 
wards dusty foreground sources, and Pchance would be 
quite inappropriate for this sample. However, given the 
highly confined distribution of dust in local galaxies and 
the observed density of galaxies on the sky, it would be 
surprising if a large fraction of GRB sightlines turned 
out to be attenu ated; indeed, more detailed analysis by 
iSudilovskv et al.l ([2009) has also recently shown that it 
cannot explain the GRB/QSO discrepancy either. 

For the purposes of this paper, we will assume no par- 
ticular bias in GRB or dark GRB sightlines. We shall 
return to this issue when discussing the implications of 
our large putative detection fraction in fJS] 

3.3. Host Photometry 

We used aperture photometry within IRAF to measure 
the flux of all candidate host galaxies, using a 1.0" aper- 
ture in all cases except for GRB 080319C, whose host 
is highly extended and a 2.0" aperture was used. In a 
few cases, the afterglow position was within the outer 
point-spread function (PSF) of a bright star, which was 
subtracted prior to photometry using various techniques 
(depending on proximity and brightness, discussed be- 
low) to avoid the complication of a variable sky back- 
ground as discussed in the relevant sections below. The 
resulting aperture magnitudes are presented in Table [3l 
A false-color mosaic of all imaging observations is pre- 
sented in Figure [TJ 

3.4. Infrared Observations 

Two events in the sample, GRBs 061222 and 070521A, 
are of particular interest. Both events were extremely 
X-ray bright, were not detected optically, and were 
observed at infrared wavelengths with large telescopes 
within a few hours after the burst. 

GRB 061222A was observed (ICenko fc Foxll2006H l to 
have a faint, fading IR afterglow. We returned to this 
field on 2009-05-31 with NIRC on Keck I and integrated 
for 10 exposures of 100s each in H- and iiT-bands. (5 sec 
X 20 coadds). Images were processed and stacked us- 
ing a modified Python/pyraf script originally written by 
D. Kaplan and aligned to our LRIS imaging. The field 
was calibrated using a single Two-Micron All Sky Sur- 
vey (2MASS) star within the NIRC field of view (2MASS 
J23530271-I-4632187). We detect a faint source near the 
detection limit close to but not coincident with the in- 
frared afterglow (likely a foreground galaxy very near the 
line of sight: see l4.10p . No source coincident with the IR 
transient is detected. Measurements and limits are re- 
ported in Table [3l 

GRB 070521 was observed less than two hours after the 
burst by NIRI on Gemini-North (|Cenko. Price, fc Bergeil 
|2007| ) and the lack of an IR detection imposes the deepest 
limit on a counterpart of any event in our sample. The fi- 
nal UVOT-calibrated XRT position contains a red source 
(well-detected in K and H, weakly detected in RG850, I 
and V) near the eastern edge. To rule out variability of 
this source, we acquired 24 x 60 s exposures in i^-band 
and 18 x 60 s exposures in iJ-band on Gemini-North on 
2009-02-01 (UT), 2.5 years after the burst. The object is 
still detected in this imaging with no evidence for fading 
photometrically or in image subtraction of the frames. 



TABLE 3 
Keck Host Galaxy Observations 



GRB Field Obj. 


-^chance 


Filter 


Magnitude^ 


AB magnitude' 


c\'^f^A.^9 A 


0.06 


9 




94 04 + n 04 






R 


99 1 ZL -1- n 09 


99 9fi + n 09 




0.45 


9 
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R 
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\^ 




9 
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R 








0.40 


9 


> 27.05 


> 26.98 






R 


94 «c: _l_ n 94 


94 07 + n 24 


^'^f^d^ r A 

UvJU'iJ-UiT- 


0.005 


9 


9/1 1 1 _|_ n 0*^ 


94 on 0*^ 

Z'i.UU ^ U.UO 






-fx 
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9*3 1 n _L n ri9 
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9 


> 25.0 
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R 


> 24.8 
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n^riTi 3 A 
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9 
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R 
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9 
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R 
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9 
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R 


9Q 4fi _|_ n OA 
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c 


0.22 


9 
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24 ^4 n 04 






R 


9Q Q7 _|_ n O'^ 


94 OS n 0^^ 

Z'i.UO ^ U.UO 




u.uu 




\ 97 9 


\ 9R S9 






y 
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J 


94 _|_ n 94 
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J 
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zo.zu ^ u.oo 
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zo. -LU ^ u.oo 






}j 
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K 


> 22.23 


> 24 03 




0.02 
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25.25 ±0.17 


25.20 ±0.17 






z 


24.10 ±0.16 


24.04 ±0.16 






J 


22.52 ±0.20 


23.40 ±0.20 






H 


21.58 ±0.09 


22.94 ±0.09 






K 


20.95 ±0.10 


22.78 ± 0.10 


080319A 


0.03 


9 


24.63 ±0.03 


24.58 ± 0.03 






R 


23.85 ±0.06 


23.98 ± 0.06 


080319C 


0.01 


9 


23.08 ±0.03 


22.99 ± 0.03 






R 


22.22 ±0.03 


22.32 ± 0.03 


080320 




9 


> 27.25 


> 27.20 






I 


> 25.3 


> 25.70 



Not corrected for Galactic extinction. 
^ Corrected for Galactic extinction. 

The final IR photometry is presented alongside the opti- 
cal photometry in Tabled 

3.5. Spectroscopy 

In several cases bright host candidates without af- 
terglow absorption redshifts available were suitable for 
spectroscopic follow-up. All spectroscopic integrations 
were conducted with longslit spectroscopy on LRIS, us- 
ing the 400/8500 grating (red side) and 600/4000 grism 
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08031 9A 



080319C 



080320 



Fig. 1. — False-color mosaic of all 14 dark GRB host fields using our Keck LRIS imaging acquired between 2005 and 2009. The 90% 
confidence afterglow positions overplotted in each case. X-Ray (XRT) error circles are cyan-colored, optical positions are green, and infrared 
positions are red. All images are 11.8" on each side with north towards the top and east to the left. See Table[2]for filter information. In 
most cases images are constructed using two filters, with the green channel interpolated using a geometric mean. 



(blue side) with the D560 dichroic, giving continuous 
spectroscopic coverage from the atmospheric cutoff to 
9200A(using the old LRIS red chip) or out to 10400A (us- 
ing the new LRIS red chip, which has greater quantum 
efficiency beyond 9000 Aand improved spectral range). 
The exposures were reduced in IRAF using standard 
techniques and flux-calibrated using observations of stan- 
dard stars BD-f 262606 and BD-H174708 (red side) and 
BD-l-284211 (blue side) at similar airmass. Absolute flux 
scales were then derived using the photometry derived 
from our previous imaging. A summary of these obser- 
vations is presented in Table H) 



TABLE 4 

Keck Spectroscopic Observations 



Field 


Obs. Date 


Exp. 


Air- 


Slit 


PA 


A 




(UT) 




mass 


(") 


(deg) 


(A) 


050412 


2007-12-13 


2x900 


1.16 


1.0 


142.85 


3500-9150 


060805A 


2009-06-25 


2x900 


1.21 


0.7 


30.40 


3500-10400 


061222A 


2007-10-09 


2x1800 


1.12 


1.0 


142.19 


3500-9350 


080319A 


2009-06-25 


2x900 


1.31 


0.7 


105.10 


3500-10400 
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Fig. 2.— Keck/LRIS R-band image of the vicinity of GRB 
050412 showing the four host galaxy candidates near the edge of 
the error circle. The XRT error circle is relatively large; only object 
A is a statistically significant association (Pchancc < 0.1). 

3.6. Photometric redshift limits 

Even in the absence of spectroscopy, it is possible to 
place limiting a redshift on host galaxy candidates using 
the color observed our optical imaging. Absorption of 
host-galaxy continuum light from hydrogen gas in the 
ISM at either the Lyman break (912 A) or Lyman-a 
(1216 A;[G unn fc Petersonl[l965f ) will strongly suppress 
the observed flux once these features enter the g-band 
at about z = 3.4 and z = 2.3, respectively, greatly red- 
dening the g — R color and allowing us to translate an 
observed color into a limiting redshift. We assume a 
strongly star -forming galaxy templ ate (the Irr template 
from hyperz [Bolzonella et al.ll2000l | , which due to its in- 
trinsic blueness provides the most conservative choice) 
with no internal extinction, then apply a simple IGM at- 
tenuation correction from .Madau (,1995. ) to measure how 
its observed g — R color evolves with redshift. At suf- 
ficiently high redshifts, the Lyman-a forest and Lyman 
break sufhciently redden the galaxy light enough to be 
inconsistent with observations, generating a simple limit- 
ing photo-z. If the redshift is known or well-constrained, 
a similar procedure can also be used to limit the internal 
extinction Ay. 

4. DARK BURSTS AND HOST GALAXIES 

4.1. GRB 050412 

The gamma-ray light curve of GRB 050412 shows no 
unusual features with a single peak and a long tail, and 
the prom pt emission fluence (15-200 keV) is 9.6 x 10~^ 
erg cm" -2 (jTueller et al.ll2005l ). near the median value for 
S'zwj/t bursts. The X-ray counterpart, however, is highly 
unusual. Swift slewed to the position after only 99 s, and 
detected a fading source inside the BAT location. How- 
ever, the X-ray flux (after decaying slowly in the first 100 
s of the exposure) plummeted abruptly starting around 



300 s, with a decay index (defined by F (x t ") of a ~ 3, 
and w as not detected after about 1200 s (jMineo et al.l 
[2001 . A Chandra X -ray Observatory Target of Oppor- 
tunity observation at 5 d (jBerger fc Foxl l2005f) failed to 
detect the counterpart. 

This burst was relatively well-positioned for ground- 
based follow-up, and was tracked by several telescopes 
including P60, all of which failed to identify a fading 
counterpart. Two additional observations deserve partic- 
ular note: a Subaru integration at 2.3 hr whic h identified 
no afterglow to i? > 24.9 mag (iKosugi et al.| [2005). and 
rapid PAIRITEL follow-up which identified no infrared 
afterglow in observations starting at 175 seconds after 
the burst trigger. Nondetections at such early times are 
rare among PAIRITEL- followed bursts (B. Cobb et al., 
i n preparation ) . 

iMineo et al.l (|2007t ) speculate that the lack of after- 
glow flux of GRB 050412 might be the result of an ex- 
tremely low-density environment suppressing the after- 
glow flux: a "naked" burst. In this case, the X-ray after- 
glow is interpreted as being completely absent, with the 
sharply-decaying light curve attributed to photons from 
high latitude from the burst itself whose arrival at Earth 
is del ayed by the curvature eS'ect (|Kumar fc Panaitescul 
|2000() . A handful of other si milar even ts exist in the 
literature, as discussed by Vetere et al.l (2008). How- 
ever, such events are very rare (at most a few per- 
cent of burst s): plotting gamma-ray fluence v ersus 
X-ray flux ( Gehrels^eLaLl 120081 : iNvsewander et al.l 120091 : 
iPerlev et al.li2009ai ). 050412 is one of only a handful of 
outliers with extremely low X-ray to gamma-ray ratios. 
The optical and IR nondetections are quite consistent 
with this picture — indeed, in terms of /3ox, the available 
constraint of /? < 0.5 — 1.0 is nothing unusual. The dark- 
ness appears to be intrinsic, not due to absorption. 

Presumably because of the weak and short-lived X-ray 
detection, the error circle of this event is large. A UVOT- 
corrected XRT position is not avail able, so the best avai l- 
able position is the one reported bv lMoretti et al.l (120061 ): 
a = 12:04:25.19, S = -01:12:00.4 (unc. 4.2")^^. 

A total of four sources are located within this region in 
our imaging, all of which are on the edge of the error cir- 
cle (Figure [2]). The first object (A), which was reported 
by several groups in the GCN circulars (jJensen et al.l 
[20051: iFox et al.llMs '). is bright [Rab = 22.3 ± 0.17)2° 
and very red {g — Rab ~ 1-8) with no clear emission 
lines over our spectral range in spite of its continuum 
brightness, which may suggest that it is an old galaxy 
with little star formation at moderate redshift (alterna- 
tively, it may also be an extremely luminous galaxy at 
2.3 > z > 1.4). Fitting line templates to the spectrum 
results in a best-fit redshift of z ~ 0.6, but this is based 
on low-S/N absorption features. In spite of the large 
XRT error circle, the brightness of the source gives a low 
^chance of 0.06, making this a probable (though by no 
means definitive) host candidate. 

Several additional, much fainter objects are also 
present near the edge of the XRT error circle. One neigh- 

All positional uncertainties in this paper are reported as 90% 
confidence error circles. 

All reported host AB magnitudes and colors are corrected for 
Galactic extinction. Afterglow magnitudes or those quoted from 
other sources are in the original reference system (Vega if BVRI, 
SDSS if griz) and are not corrected for extinction. 
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boring source (B) is not reported in any circular (likely 
because it was outside the original XRT error circle in 
the GCN circulars). It is marginally detected in both 
filters (i?AB = 25.2 ± 0.4) and has a rather typical color. 
A third source (C) was noted in the Subaru imaging of 
iKosugi et al.l ()2005l ) as being near the center of the orig- 
inal GCN XRT error circle. It is weakly detected in our 
g-band imaging {g = 25.84 ± 0.18) and marginally de- 
tected in our i?-band imaging {Rab — 25.46±0.39) which 
is consistent with the report of a marginal detection with 
R w 26.0 in the Subaru imaging. Finally, a fourth source 
(D) is at the top of the error circle and is detected with 
significance greater than 2a in R-band only. It is very 
red, with g — Rab > 2 mag. All three of these additional 



sources have R 



chance 



values of order unity. 



The large XRT error circle, and the fact that all avail- 
able host candidates are near its edge, makes host assign- 
ment particularly difficult in this case. The only object 
whose presence in or near the error circle cannot be at- 
tributed to a chance alignment with probability of order 
unity is the brightest one (source A), but especially given 
that the original XRT position did not even include this 
source there is plenty of reason to be skeptical about the 
association. If this is indeed the associated object, the 
combination of its red color, lack of lines, and perhaps 
even the fact that it is nearly outside the XRT error 
circle is particularly intriguing given the possibility of a 
very low circumburst density indicated by the X-ray light 
curve. 

4.2. GRB 050416A 

GRB 050416A (actually an XRF) is the second-lowest- 
redshift event in the P60 sample. This GRB did have 
an optical afterglow that was detected by P60 and many 
other telescopes — including the UVOT in its ultraviolet 
filters, suggesting that while this is a dark burst, it is 
perhaps a borderline case. Indeed, in terms of f3ox (equal 
to 0.37 for this burst) this event is only slightly under the 
Jakobsson criterion. 

The afterglow of GRB 05041 6A has been stud- 
ied in detail by many authors (iHolland et al.l 120071 : 
iMangano et"aH l2007t ISoderberg et al.l l2007| ) "and the 
presence of line-of-sight dust which may contribute to 
its optical faint n ess is , in principle, well-constrained. 
ISoderberg et al.l (|2007f l estimate Ay ~ 0.87 (using a 
Milky Way template), which compared t o the majority 
of GR Bs is already quite high, although IHolland et al.l 
()2007f ) derive a significantly lower value of Ay = 0. 24. 
This lower value is also favored by iKann et "all ()2007f ) . 

The host galaxy color is moderately red: q—Rab = 0.8; 
in part this is likely due to the presence of the 4000 A 
break between the q and R band s at th e emission redshift 
of z = 0.6535. ISoderberg et all ()2007f ) detected the host 
in the HST F775W filter and estimate / = 22.7 ± 0.1, 
corresponding to a significantly bluer color of (i?— /)ab ~ 
0.15. Neither of these values constrain the host extinc- 
tion strongly. However, on the basis of the observed emis- 
sion line ratio of H^/ Hp = 0.3 ±0.1, they conclude that 
the host galaxy does likely harbor significant extinction. 

4.3. GRB 050607 

GRB 050607 is at the faint e nd of Swift GR Bs, with a 
fluence of 8.9 x 10""^ erg cm-^ (jRetter et al.|[2005l Un- 
fortunately, optical follow-up of this burst was greatly 



complicated by the presence of a bright {R « 16) star 
only 4" away from the burst location. As a result, the 
P60 imaging of this burst is quite shallow, and no after- 
glow was detected in any filter. However, even if stellar 
contamination were not a problem it is unlikely that P60 
would have detected the afterglow, since m uch deeper ob- 
serva tions with the KPNO 4m telescope (jRhoads et al.l 
|2005() do detect a transient with / = 21.5 at 10 minutes, 
below the typ i cal P6 limit even in an uncrowded field. 
iRhoads et al.l ()2005l ) also note that the optical color is 
quite red, with /3opt > 1-5: suggesting either substantial 
dust extinction or a high redshift {z = 3-4). 

The bright nearby star that complicated the P60 fol- 
lowup causes substantial difficulties for host follow-up 
also. The star is saturated in our imaging, making PSF 
subtraction difficult, and the crowded field leaves no 
bright isolated template stars with which to accurately 
measure the PSF. We fit and subtract the PSF of the 
nearby star (excl uding the saturated core) using galfit 
(iPeng et al.ll2002l ). and identify no obvious source at the 
position of the optical transient. Therefore we are un- 
able to strongly distinguish between the extinction and 
high-redshift possibilities, though the _B-band afterglow 
detection imposes a limit of about z < 4. 

4.4. GRB 050713A 

GRB 050713A is another well-st udied burst— n iainly 
at X-ray and h i gher energies ([Morris "eTall I2OO7I : 
iGuetta et aP 120071 : 1 Albert et al.ll2006l ). as unfortunatelv 
the optical coverage is much more limited. It is bright, 
near the top end of the Swift sample in both gamma-ray 
and X-ray flux. The associated optical afterglow, how- 
ever, is quite faint: RAPTOR triggered on this burst 
and observed the event towards the end of the gamma- 
ray emission, but even at that point the event was only 
marginally detected with a peak magnitude of i? « 18.4 
(Wren et_al. 2005). Several prompt-emission flares at 
this time are seen in the X-ray and not the optical, but 
even after the X-ray flaring subsides the optical-to-X- 
ray index remains shallow at /3ox ^ 0.3. Unfortunately, 
this afterglow was detected in only R and / filters^^ and 
as a result the optical slope is only poorly constrained 
(/3opt — 1-4 ± 1.0) and on its own docs not constrain the 
redshift of or extinction towards this GRB. 

The position of this GRB is within the outer halo of 
a extremely bright star (1.1' from HD 204408, V 6.6 
mag). As a result, the region of the GRB is mildly com- 
promised by a variable background, which we remove by 
applying a median filter over the region of the image 
around the GRB position. After this step a source coin- 
cident with the optical position is clearly visible in R and 
marginally detected in g. The color of g — -Rab — 0.4 ± 
0.3 does not constrain the nature of the galaxy given the 
unknown redshift. It docs limit the redshift to z < 3.6, 
ruling out any contribution of Lyman absorption to the 
observed afterglow faintness. 



4.5. GRB 050915A 

GRB 050915A is genuinely dark by all definitions. It 
was followed up rapidly by several instruments, but only 



•^^ Detections in JHK filters have been reported bv lHeartv et al.l 
120051 but the photometry has not been made pubhc. 
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detected by one: the robotic infrared telescope PAIRI- 
TEL, which marginally detected a transient in iJ-band 
only {H = 18.25 ± 0.16). This is contemporaneous with 
R and /-band nondetections with the P60 that require 
an afterglow spectral index of about Popt > 1.45, outside 
the range observed for typical unextinguished afterglows 
but only weakly constraining on the rest-frame extinc- 
tion without additional constraints on the redshift and 
spectral index. Furthermore, although this is not a par- 
ticularly bright event in X-rays or gamma-rays, (3ox is 
clearly below the canonical dark value of 0.5. There is 
no evidence of X-ray flaring or a flat energy reinjection 
phase after about 100s. 

A faint ga laxy, previously discovered by 
lOvaldsen et al.l (|2007), is well-detected consistent 
with the XRT position in all filters in which it was 
observed {g, V, R, and I). It is somewhat offset (by 
about 1.1") from the IR position, although because of 
the relatively low-significance detection of the infrared 
afterglow the 90% confidence circle is large and its edge 
skirts that of the optical disk. While Pchanco is still low 
(0.06), we admit that this is one of the more tenuous 
associations in the sample. 

While the optical detection of the host alone rules out 
a high-redshift origin, VLT spectroscopy of this galaxy 
(P. Jakobsson et al. in preparation) has revealed a sur- 
prisingly low redshift of z ~ 0.4, indicating an extremely 
underluminous system (Mv(yiB) « —17.4) and requiring 
a significant (though not, in this case, particularly large) 
dust column to explain the redness of optical afterglow. 
Consistency of the combined X-ray and optical data re- 
quires Ay ^0.5 mag independent of extinction law. 

The blue colors of this galaxy indicate a young popula- 
tion free of widespread dust (global Ay ^1-0 mag from 
our template modeling). This limit is not inconsistent 
with the relatively modest minimum extinction inferred 
from the afterglow. 

4.6. GRB 060210 

GRB 060210 provides significant insight into the 
dark burst phenomenon. The optical afterglow of this 
burst was fairly bright, but only in the reddest bands 
(i? and /), peaking around 19.5 mag at a relatively 
late time of 600s following an extended episode of X- 
ray and optical flaring. Afterglow spectroscopy by 
ICucchiara. Fox, fc Bergerl ((2006) confirmed that this is a 
(moderately) high-redshift event at z = 3.91, explaining 
the steep fall-off towards the optical bands. In addition, 
there is significant evidence for high-redshift dust, given 
that even optical filt ers redward of Lyman-g are sig- 
nificantly suppressed (jCurran et al.ll2007f l. ICenko et al.l 
(|2009D estimate A y = 1.21^°;^^ mag (in agreement with 
iKann et al.l l2007l ). which at the burst redshift corre- 
sponds to ~ 4 mag of extinction in the observed i?-band 
using an SMC template. 

We imaged the field on two occasions; a relatively short 
R integration followed by deeper g and / observations. 
Unsurprisingly, nothing is detected in g-hand, which falls 
below the wavelength of Lyman-a and is likely to be 
heavily obscured. However, a bright source is detected 
at the OT position in R and I.^^ The offset between 
this object and the OT is less than 0.5" (Pchancc < 0.01) 

This is not the object mentioned in lHeartv et al.ll2006l . which 



and the association is further bolstered by the g-band 
nondetection. This therefore likely represents among the 
highest-redshift GRB host galaxies detected to date, as 
well as among the most luminous {Mn < —20.2 for a 
starburst template). In spite of the optical extinction, 
redward of Lyman-a the color of the object is quite blue, 
with {R — T)ab = O.libO.3 (the large uncertainty is dom- 
inated by the poor calibration of this field using USNO 
standards). Given that the R and / bands correspond 
to wavelengths well into the ultraviolet at this redshift 
(1300-1700 A) where dust absorption is extremely effi- 
cient, this suggests that the average observed extinction 
cannot be high, though given the lack of knowledge about 
the extinction law it is difficult to constrain this formally. 
For an assumed SMC-like extinction law, the host extinc- 
tion is Ay — 0.25±0.25, which is certainly much less than 
the inferred extinction from the afterglow. 

4.7. GRB 0605 lOB 

The spectroscopic re dshift of this event {z = 4.941, 
lPricell2006HPrice et al"]l2007>) is the highest in the sample 
and among the highest for any burst to date. At this 
redshift the Lyman-a transition is shifted well into the 
optical band, and consistent with this the flux in the 
P60 R and i bands is strongly suppressed. Blueward of 
i?-band the OT is not detected. Unfortunately, this is 
one of the few bursts which displays clear flaring activity 
in the X-ray band as late as 1000 seconds after the GRB, 
making a consistent estimate of /3ox difficult, though as 
measured in i?-band the burst is clearly dark for almost 
any assumption of the X-ray afterglow behavior. 

Optically, coverage of this burst was sparse, and both 
R and i filters are affected by Lyman-a absorption, mak- 
ing it difficult to estimate the extinction. However, 
the z — Jab color of 0.0 ± 0.4 ( based on the J-band 
point of lPrice. Cenko. fc Foxll2006f) requires Ay < 0.5 for 
Popt > and SMC-like extinction. In addition, the late- 
time Pm-x (using the J-band point) is actually 1.0 
and entirely normal, giving further evidence that the ex- 
tinction is negligible. Because of the known high redshift, 
our integration on this source was particularly long (ap- 
proximately one hour), though the quality of the images 
is poor due to bad seeing (1.4"). No object was detected 
at the P60 position or anywhere inside of XRT and XMM 
X-ray error circles in cither the Roy g filters to 26th mag- 
nitude. 

The host galaxy of this burst was imaged b y the Spitzer 
Space Telescope in a study conducted by ICharv et al.l 
(|2007i ). and successfully detected with a flux level of 
0.23 ± 0.04 yuJy. Our g-band nondetection can be in- 
terpreted as support of this association (a detection of 
a galaxy blueward of the expected Lyman break in or 
near the optical position would indicate that the Spitzer 
source was actually an intervening source at lower red- 
shift). Given the high redshift, the R nondetection is not 
surprising either; our limit of i? > 26 corresponds to a 
luminosity of Mn > —20.5, which is still consistent with 
the luminosities of the majo rity of GRB hos ts which have 
been observed to date (Fruc hter et al.l[2006l ) and with the 
sub-L* nature of the reported Spitzer host (jCharv et al.l 
I2OOI . 

according to that note is 2-3"north of the XRT position. No source 
is detected at that position in our imaging. 
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Fig. 3.— R-band imaging of the field around GRB 060805A. Two 
galaxies are located inside the XRT error circle with additional 
objects nearby. 



4.8. GRB 060805A 

GRB 060805 A was an extremely faint Swift burst, with 
a fluence value in the bottom 3 percent of all Swift long 
GRBs (the burst was not detected at all above 100 keV). 
The X-ray afterglow is extremely faint: « 3 x lO^'' mJy 
even at 100 seconds. 

From this perspective it is no surprise that P60 (and 
all other optical instruments) failed to detect an opti- 
cal afterglow, and indeed the limit on the optical to X- 
ray slope is effectively nonconstraining at /3ox < 0.7. 
The low observed flux and fluence suggest an intrinsi- 
cally low-luminosity event, though a typical- luminosity 
GRB at sufficiently high rcdshift could also appear faint 
simply because of its great distance. Our imaging obser- 
vations favor the former interpretation: two host galaxy 
candidates are present within the XRT error circle: one 
bright object (object "A" of Figure [3l Rab ^ 23.6) at 
the southwestern edge and a second, fainter source (ob- 
ject "B", i?AB ~ 24.6) slightly northeast of center. The 
colors are significantly different: the brighter source is 
blue with g — Rab ~ 0; the fainter one is redder with 
g — Rab ~ 0.8. Unfortunately, we are not able to distin- 
guish which is the correct host given the size of the XRT 
error circle. 

Our spectroscopic observation of this source used a slit 
angle covering both sources (A and B) . Only the brighter 
object (A) shows a noticeable continuum trace in our 2D 
spectra. No line features are observed over the spectral 
range down to the atmospheric cutoff; the nondetection 
of Lyman alpha or associated absorption features implies 
approximately z < 1.8. The nondetection of Lyman al- 
pha at the position of object A may impose a similar 
redshift constraint on this object also, but this conclu- 
sion is less robust. The redshift limit implied by the g — R 
color is z < 3.8. 
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Fig. 4. — Stacked V- and 7-band image of the field near GRB 
060923A. A faint galaxy is marginally detected coinci dent with the 
bright est, central region of the galaxy as also noted bv lTanvir et al.l 
II2008I ). A projection from the galaxy appears to extend towards 
the southwest. 

4.9. GRB 060923A 

One of the clearest examples of a dark burst in the 
sample is GRB 060923A. Though not a particularly high- 
fluence event in gamma-rays or in X-rays, this burst was 
observed very early in the NIR (< 1 hr) using UKIRT 
(jTanvir et al.l l2006f) and shortly thereafter with both 



Keck and Gemini (jFox. Rau. fc Of ek 2006) . A transient 
was detected in ii'-band in all of these observations, but 
not in any bluer filter including J or H . One possible ex- 
planation for this would be an extremely luminous event 
at high rcdshift [z > 15). However, later optical follow- 
up by Tanvir et al. (2008) identified a host galaxy exactly 
coincident with the IR location, marginally detected in 
our imaging as well in V- and /-bands (Fig ure ]2) ■ It is 
not detected in B or i?G850. iTanvir et aLl (j2008) esti- 
mate that for z = 2.8 about Ay ~ 2.6 would be sufficient 
to explain the inferred absorption. 

The host galaxy is fairly but not remarkably red in 
the observed-frame optical: {B — V)ab ^ 0.5 and 
{V ~ I)ab = 1.0 ± 0.4. A nondetection in RG850 rules 
out continuation of this trend further to the red, im- 
plying that the spectral energy distribution (SED) fiat- 
tens towards the rest-frame o ptical, incons i stent with a 
highly dust-obscured source. [Tanvir et all (|2008f) addi- 
tionally report {R — K)ab ^ 2.1, which is not unusual for 
moderate-redshift galaxies. We attempted to fit model 
SEDs using the combined BVRIzK photometry, but due 
to the poor detections in all filters no reliable model con- 
verged. Further, only a redshift l imit of .z < 4.4 is p ossi- 
ble from our photometry, though [Tanvir et ani2008l con- 
clude that z < 4.0 based on the combined properties of 
the X-ray and optical afterglows.] Additional photome- 
try will be necessary to reliably constrain the extent of 
extinction and other properties of the host, but as with 
most other galaxies in our sample the host-galaxy pho- 
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Fig. 5.— y-band image of the field near GRB 061222A. Two 
objects with similar magnitudes and colors are slightly blended; 
only the northern o bject (A) is consisten t with the Gemini infrared 
afterglow position llCenko &: Fox|[2006al ). The position of the slit 
used in acquiring spectroscopy of the two sources is also shown. 



tometry does not demand large amounts of dust. 

4.10. GRB 06 1222 A 

At high energies, GRB 061222A is among the bright- 
est events in the sample. The gamma-ray light curve 
contains numerous separate pulses and extensive flaring 
out to '^100 s, and the X-ray flux is also bright, well- 
detected by the XRT out to 10^ s. As measured at ^ 11 
hr the X-ray flux from this event is in the top 2% of all 
Swift GRBs. 

Several other telescopes in addition to the P60 ob- 
served this event at early times, generally obtaining rel- 
atively shallow limits. However, NIRI was triggered at 
Gemini in K-h&nd only (jCenko &: Foxl2006af) . and a faint 
source was identified that later faded, confirmin g this to 
be an infrared afterglow (|Cenko k, Foxl [2006bl ) . Unfor- 
tunately no deep imaging was acquired in other filters. 
Howe ver, this event was also detected in radio using the 
VLA (jGhandra fc FrailllMol . 

Two blended but seemingly distinct sources are ob- 
served near the afterglow position (Figure[5]): one (source 

A) coincident with the IR transient and a second (source 

B) offset by about 1" to the southeast. We identify the 
former as the host galaxy. The two objects have similar 
colors, though photometry is complicated by the close 
blending, especially in the redder filters where neither 
object is well-detected. Both galaxies are quite blue, 
with {B — F)ab 0.0 mag, {V — /)ab ~ 0.5 mag, 
(/ — 2;)ab ^ 0.3 mag. Only object B is detected in the 
infrared, but both galaxies are clearly very blue in IR 
colors as well: for object A, (/ — i^)AB < 1.0 mag; for 
object B, (/ - K)as. = 0.8 ± 0.3 mag. 

Our LRIS longslit spectroscopic observation placed 
both objects along the 1" slit for two exposures of 1800s 
each. The telescope was dithered 5" between the ex- 



posures. The blue-side exposure was reduced normally, 
though the severe fringing on the red side was only 
removed effectively by subtracting the two exposures, 
which cleanly removed the fringe pattern. We extracted 
spectra separately for both sources (A and B) along the 
slit near the afterglow position. Interestingly, despite 
similar colors these galaxies are not at the same redshift. 
The fainter, northern object (A), which we identify as 
the host galaxy, has a strong emission line at 3758 A. No 
flux is observed at this position in the southern object 
(B). At the same time, between two sky lines on the red 
side another bright emission line is clearly observed at 
8014 A in this case consistent only with the position of 
object B. The spectra and putative lines of both objects 
are shown in Figure [6l 

The strong line in the blue part of the host-galaxy 
spectrum strongly suggests Lyman-a at a redshift of 
z = 2.088. An alternate possibility is [OH] at z = 0.008, 
but this would require an extraordinarily small and un- 
derluminous galaxy as well as imply the presence of Ha 
at 6617 A, which is not observed. Galaxy B cannot be 
at this redshift — its solitary line, if interpreted as [OH], 
indicates z = 1.151. (Alternatively, the line could be as- 
sociated with Ha at z = 0.22, but this would predict the 
presence of [OH] at 4550 A which is not observed.) 

At the observed redshift, any suppression of the bright 
optical afterglow predicted by the bright X-ray counter- 
part must be due to dust extinction. The darkness of 
this burst is truly extreme: even in the observed K- 
band, approximately 4 mag of extinction are necessary 
if we assume the minimum synchrotron intrinsic spec- 
tral index of /3ox = 0.5. At the observed host-galaxy 
redshift oi z — 2.088, this corresponds to approximately 
Av > 5.0 mag (nearly independent of the choice of ex- 
tinction law).^^ 

Given the enormous amount of extinction inferred from 
the faint infrared afterglow, one might expect that the 
relative amount of extinction in the observed optical 
bandpasses should be even greater — yet the host can- 
didate is relatively bright {V ^ 24 mag) and extremely 
blue, showing no signs of reddening at all: the broadband 
color strictly limits the host-galaxy Ay < 0.5 mag. 

4.11. GRB 070521 

Like GRB 061222A, GRB 070521 was a bright GRB 
with a bright X-ray counterpart. In addition to the stan- 
dard P60 foll ow-up, observat ions commenced at P200 
within 1 hr ()Rau et all I2007D and at both Keck and 
Ge mini (including, in the latter case, JHKs IR imag- 
ing: iMinezaki fc Pric3 12007) within 2 hr. As described 
in §3.4( no transient source within or near the XRT er- 
ror circle was identified in any of this imaging despite 
the rapidity, depth, and relatively long wavelengths of 
these observations, making this burst the darkest in the 
sample. 

In our observations, the most recent UVOT-enhanced 
XRT error circle includes a red, pointlike object near its 
eastern edge (Figure [7]) . It is strongly detected in the 

■^•^ It is conceivable that the foreground object may also con- 
tribute to the extinction, but the blue colors of both this foreground 
object and the host (which would be reddened by a similar degree 
as the afterglow) make it unlikely to be a large contributing factor 
to the large absorption demanded by the afterglow. 
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Fig. 6. — LRIS red-side spectrum of the host galaxy GRB 061222A and a nearby (in projection) object placed along the slit (Figure |5j 
with insets showing detected emission lines. A strong line is observed in the host galaxy (object A; top) at 3758A which we interpret as 
Lyman-a at a redshift of 2 = 2.088. No other objects are observed over our spectral range. Despite the small offset and similar broadband 
color, object B is not at the same redshift. No flux is observed at the location of the putative Lyman-a line; instead, we detect a single 
line at 8015A which we interpret as the [OII]3727 doublet at a redshift of z = 1.151. 



NIR filters (except J, which was a relatively short expo- 
sure). However, in /-band it is only marginally detected, 
slightly blended with another source located outside the 
XRT error circle, and was only detected in our y-band 
imaging after a second visit to the field: uniquely among 
the host-galaxy candidates in this sample, this object is 
quite red. No other objects are present within the error 
circle at either optical or infrared wavelengths. 

Thanks to the large suite of broadband photometry 
available for this object, we have been able to model the 
host SED and estimate an appro ximate photometric red - 
shift. Using the package hyper z (jBolzonella et al.|[200ClD . 
the SED is well fit by a late-type galaxy template at a 
redshift of z = 1.35^q'32 with a stellar age of 360 Myr 
and an extinction of only Ay = 0.4 mag. The apparent 
redness is, therefore, more likely to be due to the pres- 
ence of the 4000 A break rather than dust: indeed, the 
JHK SED redward of this break is quite normal. There- 



fore, as with the other host galaxies in our sample, little 
dust extinction is demanded by the host data. 

The amount of extinction required by the afterglow of 
this burst is as phenomenal as for 061222A. Assuming 
an intrinsic afterglow /3ox > 0.5 at 10^ sec, the deep 
Gemini limit requires an extinction of at least 4.7 mag 
in the observed if -band. At the putative host redshift of 
z — 1.35 this corresponds to a limit of Ay > 9 mag (over 
the 95% confidence redshift range of z = 0.95 — 2.05, the 
constraint is Ay > 6 mag). A similarly large amount of 
extinction in the host SED is ruled out by our template 
modeling. 

4.12. GRB 080319A 

GRB 080319A was a relatively bright GRB, though 
both its X-ray and optical afterglows are unremarkable, 
and the observational coverage sparse — likely as a result 
of the intense focus on GRB 080319B which occurred 
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Fig. 7. — .ff-band imaging of the host galaxy of 070521 from NIRI on Gemini-North alongsi de J-band im aging from Keck. SI and S3 ar e 
nearby, unassociated objects that were proposed as possible hosts in the GCN circulars (Hatt ori. Aoki fc Kawai 200^: iPerlev et ani2007l ). 



only 27 minutes afterward in the same part of the sky. 
Swift 's initial slew to this burst was also delayed by 500 
s due to an Earth constraint (in total, the XRT observed 
for only two epochs — at ~ 1 ksec and briefly at 4 ks). 
Optically, the afterglow is detected by P60 in Riz filters 
and in a single epoch with the UVOT at approximately 
600 s. PAIRITEL also successfully detected the afterglow 
in JHK before slewing to 0803 19B. The IR color is also 
red and consistent with the optical color, for an overall 
optical-NIR spectral index of (i — 1.5. This is suggestive 
of significant extinction. 

A relatively bright galaxy is located coincident with the 
P60 optical afterglow position. As with other galaxies in 
our sample, the optical color is not unusual [g — Rab = 
0.60 ± 0.06). While this single color does not strongly 
constrain host extinction, as with other bursts the rel- 
ative brightness of the host combined with the lack of 
obvious redness does not give any reason to suspect its 
presence. Spectroscopy reveals no line features over our 
spectral range redward of the atmospheric cutoff, limit- 
ing the redshift to z < 2.2. At this redshift and assuming 
an intrinsic spectral slope (3 < 1.2, the lower limit on the 
extinction implied by the photometric SED is Ay ~ 0.25 
mag (SMC extinction). Any deviation from these as- 
sumptions (lower redshift, shallower intrinsic slope, or 
other standard extinction laws) would require additional 
extinction, implying a lower limit on the extinction of 
Ay > 0.25 mag. 

4.13. GRB 080319C 

GRB 080319C was a bright, hard burst, and trig- 
gered several satellites in addition to Swift inc l uding 
Suzaku, Konus, and Agile [M arisal di et al. I 120081 : 
iGolenetskii et al. Il2008t lOnda et al. 1 120081 ). The after- 
glow is relatively unremarkable at late times, and was 
detected by the UVOT in filters as blue as U and so 



clearly is not as "dark" as other objects in this sample 
(/?ox = 0.36). The bur st was in fact br ight enough for 
an absorption redshift (jWiersema et al.l [20081 to be ac- 
quired, placing the event at z = 1.95. However, as is 
the case with the other bursts in the sample, the ob- 
served optical fluxes are suppressed relative to the X-ray 
flux and show evidence of reddening, which can be esti- 
mated with precision thanks to the known redshift and 
larg e numbers of filter s {Ay = 0.67±0.06 mag, consistent 
with lKann et al.ll2007T ). The optical and X-ray afterglows 
both show a dramatic flare around 200 s, after which the 
afterglow appears to decay relatively uniformly, though 
coverage is sparse. 

The host galaxy of this event is remarkably bright: 
i?AB — 22.3 mag. In fact, this value is consistent with 
the reported P60 flux from observations taken the night 
after the GRB and was likely serendipitously detected 
even by this small-aperture telescope. The relatively high 
redshift makes this particularly remarkable: the absolute 
magnitude of this galaxy for a flat-sp ectrum fc-corr ection 
IS Mr = -22.6mag(~4i» at z ~ 2: lReddv et "al][2 008). 
which would make it the second most luminous GRB host 
galaxy known (s econd to the even mor e remarkable host 
of GRB 081008. ICucchiara et al.|[2008l if its reported lu- 
minosity is real.). The color is blue (5 — -Rab = 0.33±0.05 
mag). At the observed redshift this is not strongly con- 
straining on the host dust: the 2175 A bump (if present — 
there is no evidence for it in the afterglow SED, though it 
is not strongly constrained) would shift into the i?-band 
and as a result the broadband extinction is essentially 
gray around these wavelengths. 

An alternate hypothesis for the galaxy detected in our 
imaging associ ates it with the z = 0.8104 Mg II absorber 
in Table 54 of iFvnbo etlHI (|2009l l. rather than the true 
host. The apparent brightness of the object would be 
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much less remarkable at this redshift, especially if the 
flux were combined with that of a true background host 
along the line of sight. Spectroscopy or high-resolution 
imaging of this system will be needed to confirm or rule 
out this possibility. 

4.14. GRB 080320 

GRB 080320 is a relatively faint S'wi/t burst with a 
mostly featureless light curve, though the X-ray light 
curve shows significant flaring ending at around 1000 s. 
Due to the nearly full moon and the attention towards 
the previous night's GRB080319B, the optical afterglow 
was observed only sparsely. This makes it difficult to 
accurately construct an SED of this event. However, as- 
suming no dramatic color changes or late-time optical 
flaring, all data are consistent with a very red afterglow 
color: using contemporaneous or near-contemporaneous 
epochs we estimate i — r > 1.1 mag, 2; — i « 0.8 ± 0.2 
mag, and Jab — « ~ 2.2 mag. Alone, these observa- 
tions are not sufficient to distinguish between a highly 
extinguished or high-redshift counterpart, though there 
is suggestion that both probably contribute: the SED is 
red across many filters, which is characteristic of extinc- 
tion but less so of a Lyman break. However, the J-band 
is probably not strongly suppressed relative to the X-rays 
(/3ox ~ 0.5 as measured from J-band), and furthermore 
our early-time PAIRITEL limits on this event show no 
evidence for a bright JC-band afterglow that may be ex- 
pected if this redness carried into the optical. The z-band 
detection imposes an upper limit on the redshift of 2: < 7. 

Consistent with this interpretation, we do not detect 
any host galaxy at the position of the optical transient 
to deep limits. While in principle this could simply be 
the result of a low- luminosity host, the Nh column mea- 
sured by the XRT is relatively low in comparison with 
the dark bursts in our sample for which we infer large ab- 
sorption columns, offering additional support of a mod- 
erately high-redshift origin (X-ray absorption is strongly 
wavelength-dependent, with the same column absorbing 
much more efficiently at lower energies: at higher red- 
shift these lower energies are shifted out of the XRT sen- 
sitivity window and sw amped by the Ga lactic absorp- 
tion signature — see also lGrupe et al.ir2007l ) Of course, a 
small host would predict a relatively low absorption col- 
umn as well, though significant dust extinction in such 
a system would not be expected. Nevertheless, we can- 
not rule out this scenario and can formally only place an 
upper limit on the redshift. 

5. RESULTS 

5.1. Redshift limits and the implications for high-z 

GRBs 

An afterglow detection in any optical (/-band or blue- 
ward") filter immediately rules out a high-redshift ori- 
gin.^ So does an optical detection of the host (assuming 
we have a proper identification). Using the combination 
of these two factors we can place an upper limit on the 
number of bursts in our sample which could have origi- 
nated from very high redshift (z > 7). 

Some measurable flux could be detected blueward of the 
Gunn-Peterson trough from a sufficiently luminous event at 5 < 
z <7, though such an event would show a clear photometric break. 
We find no evidence for such an event in the P60 sample among 
GRBs with unknown redshifts. 



In fact, no events of our sample are consistent with 
such a high-redshift origin. If we assume our proposed 
host associations are all correct, all 29 events in the P60 
sample have either an optical transient or an optical host 
candidate, suggesting that — in spite of Swift 's sensitivity 
and customized trigger software — it detects few events 
beyond the range that has already been probed by optical 
spectroscopy. Under this assumption, all events in the 
sample are constrained to z < 7 and all but one (GRB 
080320) to z < 5. 

Because the P60 sample is uniformly drawn from all 
Swift events, we can convert this observational statement 
to a limit on the intrinsic high-redshift fraction among 
Swifthmsts. We perform a simple Monte Carlo sim- 
ulation in which 29 events are repeatedly drawn from 
a source population with the intrinsic high-z fraction 
treated as an input parameter. To convert this to a 90% 
confidence upper limit, we then vary this input parame- 
ter until zero high-z are events are drawn in exactly 10% 
of the simulated 29-event samples (for z > 7) or zero or 
one event is drawn in exactly 10% of the samples (for 
z > 5). We conclude that, if all of our supposed associa- 
tions are correct, at most 13 percent of Swift events are 
at z > 5 and at most 7 percent are at z > 7 to within 
90% confidence. (This procedure can be generalized to 
lower redshifts also with appropriate assumptions — see 
Figure H) 

These estimates have neglected the possibility that 
some of our host associations may be chance alignments 
with foreground galaxies. To take into account the pos- 
sibility of foreground galaxy contamination, we assumed 
that 10% of high-z events in our simulation would be 
wrongly associated with a low redshift host (Pchanco — 
0.1 is the largest observed in any of our possible host 
associations) and performed the simulation again, vary- 
ing the true high-z fraction until zero apparently high-z 
events are present in 10% of the samples (for z > 7), 
or zero or one apparently high-z events are present in 
10% of the samples (for z > 5). In fact, this changes our 
constraints only slightly (by about one tenth of each per- 
centage value). We therefore conclude that, within 90% 
confidence, at most 14 percent of all Swift GRBs 
are at z > 5 and at most 7 percent are at z > 7. 

As our most conservative assumption, we may choose 
to reject two host associations completely (in spite of the 
low Pchance)- Specifically, suppose we reject two of the 
six host associations for events with no optical detection 
(for events with optical detections whether or not we 
have identified the host correctly does not significantly 
impact our conclusions about the redshift distribution) — 
that is, we assume a 33% contamination rate, which is 
much higher than that anticipated by chance. In this 
case, the data are consistent at 90% confidence with up 
to 23% of GRBs at z > 5 and up to 18% at z > 7. How- 
ever, we point out that the "weaker" associations (where 
error circles and/or offsets are large: 050412, 050915A, 
060805A) are consistent with simply being underlumi- 
nous in all bands and no more likely to be at high-redshift 
than any other burst even if their claimed host galaxies 
are unassociated. The one exception, GRB 070521, has 
independent confirming evidence for a highly-absorbed, 
low-z nature in the form of a large X-ray A'h column (as 
do the statistically firmer associations of GRB 061222A 
and GRB 060923A; see also Figure E]). 
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Fig. 8. — Cumulative redshift distribution of Swift GRBs in- 
ferred from the uniform P60 sample. The two solid gray lines 
show the absolute maximum and minimum observed redshift dis- 
tribution in the entire sample (that is, assuming all GRBs with 
unknown redshift are at the maximum possible [see Table |5] or the 
minimum possible [z = 0] redshift). The thick black line is the 
distribution omitting non-dark (/3ox > 0-5) or potentially dark 
2 < 4 events with no measured redshift (these events are, as a 
population, not likely to significantly deviate from the redshift dis- 
tribution of Swift events in general) and conservatively assumes 
GRB080320 to be at z = 6. Based on this assumption, the salmon 
region then represents statistical limits on the cumulative fraction 
of Swift GRBs originating at or below a given redshift as a func- 
tion of 2 permissible to be consistent with the observed distribution 
(10-90% confidence limits). The inferred distribution is generally 
consistent with the observed distribution of spectroscopic redshifts 
for all Swift events to date, indicating that there are no strong 
redshift biases — except possibly at the highest-2 end, where the 
observed and intrinsic rate are not as well-constrained. 

The recent detection of GRB 090423 at z = 8.2 also 
allows us to place a (relatively non-constraining) lower 
limit on the high-redshift fraction. While the POO sam- 
ple in this paper was cut off at the end of March 2008, 
P60 triggered rapidly on GRB 090423 and detected no 
afterglow to limits comparable to those discussed in this 
work. GRB 090423 was the 42nd Swift GRB on which 
P60 triggered rapidly. We perform a simple Monte Carlo 
simulation in which bursts are sampled from an intrin- 
sic population with a user-specified high-z rate, which 
is varied until a high-z event occurs as or earlier than 
the 42nd event 10% of the time. Only a rate of 0.2% 
is necessary to fulfill this criterion. Therefore, the de- 
tection of GRB 090423 requires (to >90% confidence) 
only that a minimum of 0.2% of all Swift events origi- 
nate from z > 7, which is fully consistent with our max- 
imum value inferred from the sample discussed in this 
paper. Our overall constraint on the z > 7 burst fraction 
for Swift is therefore 0.2-7 percent (to within 80% confi- 
dence). This estimate is consistent with other recent ob- 
servational limits on the high-z fraction, su ch as that of 
Ruiz-Velasco et al.ll2007l ( £ 19% at z > 6) iGrupe et all 
20071 « 7%atz > 6). and lJakobsson et al.ll200,'Tl (;7-40% 
at z > 5). 

Our results strongly constrain some theoretical mod- 
els of the evolution of the GRB rate with cosmic time. 



For example, iBromm fc Loebl (|2002f) predicted that 50% 
of all GRBs and 25% of Swift GKBs originate at z > 5, 
which we rule out. It is consistent with some more re- 
cent models that predict a low hig h-z GRB rate based on 
star formation rate (SFR) models (jBromm & Locbl200i; 
iLe fc Dermeill2007[ ) ■ luminosity indicators (.LiiOOS), and 
limits on the GRB production efficiency of Population III 
stars (iBelczvnski et al.ll2007l: iNaoz fc i3romberell2007t )— 
though some of these models predict high-z fractions 
close to our maximum value, which a larger sample may 
be able to confirm or refute. 

5.2. Constraints on dust extinction 

In Table [5] and in Figure [10] we have summarized the 
extinction constraints derived based on our host-galaxy 
redshift constraints and the properties of the afterglow. 
For most bursts in the full P60 sample, there is little 
extinction: the median Ay is about 0.5, within the range 
of values typically seen i n previous studies of optically 
well-studied bursts (e.g.. iKann et al.]l2006f ). However, 
large extinction columns are common: six bursts (out of 
22 in which useful constraints can be derived) have Ay > 
0.8 mag and three have Ay ^ 2.5 mag. In comparison, 
only two events have i?-band fluxes that are suppressed 
by hydrogen absorption at high redshift. Thanks to the 
uniform nature of this sample, we therefore are able to 
assert with reasonable confidence that the predominant 
cause of the dark burst phenomenon is dust extinction. 
Even an extinction of Ay ^ 1 mag translates to large 
i?-band extinction values at typical ^'wii/i redshifts (> 3 
mag at redshifts of z > 2). 

Unfortunately, the nature of this dust remains a mys- 
tery. The hosts of highly extinguished events tend to be 
unremarkable objects — often optically bright and with 
no evidence for large amounts of intrinsic reddening, and 
in a few cases with blue colors that appear to directly 
contradict the expectation of extremely red objects. In 
no case are the optical colors indicative of a ULIRG- 
like highly e xtinguished object, which some theoretical 
models (e.g.. lRamirez-Ruiz et al.ll2002 f) predict should be 
common among the GRB host population. This result is 
not peculiar to our study: other dark burst hosts have, 
in the large majority of cases which have been studied 
to date, also been relatively blue o bjects without clear 
photometric evidence for extinction (Jaunsen et al. ' 20081 
Ro l et al. 2007 ; Ca stro- Tirado et al. 2007; Pellizza et jj] 
I2OO6I: fCorosabel et all l2003l : iDiorgovski et all 12001 7 
though a few counterexam ples of very r ed hosts do exist 
as well (jLevan et al.ll2006HBerger et al.ll2007[ ). 

The results can be interpreted in several ways. One 
possibility is that these apparently blue galaxies are con- 
cealing their true natures: if the distribution of dust is 
sufficiently patchy, it is conceivable that what looks like 
a normal object in the optical and NIR bands could har- 
bor a massive starburst obscured from view by the same 
dust concealing the afterglow, allowing the (blue) emis- 
sion from the optically thin regions to dominate the SED 
even if they contribute little to the total SFR. Alterna- 
tively, there could be relatively little dust in the galaxy 
overall, but the GRB itself could be located deep within 
a relatively small dusty patch, such as a young molecu- 
lar cloud, though this region would have to be sufficiently 
larg e to escape the destructiv e influence of the burst it- 
self (|mxm^^Draini|2003) . A third, more exotic pos- 
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TABLE 5 

Redshift and extinction constraints on P60 GRBs 



GRB 


/3ox'' 


ijiuefeL 
AG det. 


host det. 


z 


(mag) 


ivjj excess 
(z = 0) 
(10^" cm ^) 


















< 0.49 




y ■ 


<4.5? 




<93.2 




0.35 


UVW2 


9 


6535^= 


0.87 


24 n+6 


050607 


~ 0.33 


B 


none 


<4 




<15.0 


05071 SA 


0.31 


R 


9 


<3.6 




<28.3 


05091 5A 


< 0.44 




9 


^0.4 


> 0.5 


<14.4 


060805A 


< 0.76 




9 


<3.8 




<38.0 


060210 


0.37 


R 


R 


3.91 


1.21 


8.7-1-2.1 


060q2SA 


< 0.24 


K 


Y 


<4 


~2.5 


22.1±9.2 


061222A 


< —0.19 


K 


B 


2.088° 


> 5.0 


34.6±2.8 


060510B 


0.04 


R 


3.6// 


4.941 


< 0.5 


<14.6 


070521 


< -0.10 


none 


V 


~ 1.35 


> 6 


44.1±12.7 


080319A 


0.41 


V 


g 


<2.2 


> 0.25 


<17.3 


080319C 


0.36 


u 


a 


1.95 


0.67 


<32.6 


080320 


< 0.31 


i 


none 


<7 




8.7± 3.3 


Other GRBs 














050820A 




UVWl 


n 
y 


2.615 


< 0.1 


3.4±1.5 


050908 




V 




3.35 


< 0.55 


<19.3 


060110 




R 




<5 


< 0.3 




060502A 




B 




1.51 


0.53 


<5.5 


060906 




R 




3.685 


0.2 


<31.2 


060908 




UVWl 


V 


1.884'' 


< 0.1 


<12.6 


070208 




R 




1.165 


0.96 


<38.8 


070419A 




9 


r7S 


0.97 


0.70 


<35.8 


071003 




U 




1.60435 


< 0.26 


<13.9 


071010A 




a 




0.98 


0.60 


<37.0 


071011 




V 




<5 




<60.7 


071020 




R 




2.145 


< 0.35 


<16.1 


071122 




white 




1.14 


0.58 


<10.6 


080310 




UVWl 




2.43 


0.10 


<7.9 


080319B 




UVW2 


9 


0.937 


0.07 


4.4± 2.2 


090423' 


< 0.5 


J 


none 


8.3 


0.1 


<10.6 



Note. — Properties of the full P60 sample (including non-dark bursts), modified from 
ICenko et al.l 1 120091 ) to include additional constraints based on the host galaxies and some 
additional afterglow data. The redshift can be constrained in almost all cases. See text 
for additional information. 

Only listed for bursts identified as "dark" in the sample. 

^ Extinction in the host-frame V-band along the line of sight inferred from the afterglow, 

generally assuming SMC extinct i on an d /3opt ~ 0.6. 

" Ay refere i ices: [Cenko et al.l l[2009l ): ISoderberg et al.l (|2003); ICovino et all pOOSh : 
IBloom et al] l l2009t) 

Equivalent 2 = hydrogen column density in ex cess of the Galactic value inferred 
from the X-ray spectrum, fit using the procedure of IButler fc Kocevskil Il2007al) . Only 
detections of > 2rj excess are shown; other events are displayed as upper limits (see 
Figure |9] for a less conservative assessment of Nh columns for bursts in the sample.) 
All objects for which we infer large amounts of dust extinction in the optical band also 
have unambiguous detection of excess X-ray absorption columns; no event with low or 
negligible dust extinction shows this signature. Since X-ray absorption is more efficient 
at low redshifts, this offers additional support to our association of these objects with 
moderate-2 hosts. 

Emission-line redshift. 
^ While not formally in our sample, the recent GRB 090423 is presented for reference as 
an example of a confirmed z > 7 event. Notably, this event has no host galaxy to z > 26 
a nd no sign i ficant excess A^h column uTanvir et al.ll2009ll . 

s |Dai et al.l |[20M . 

h Revised redshift from lFvnbo et al.ll2009l . 



sibility is that our templates for modeling high-redshift 
dust are incorrect, and high-redshift GRB hosts are dom- 
inated by grey extinction la ws that redden the ir stellar 
popu l ations relatively li ttle (IChen e t al. 200^ iLi et al.l 
I2008t iPerlev et al.|[200l . Unfortunately, the available 
data do not allow us to distinguish between these possi- 
bilities. 

In any case, however, our results suggest that a sig- 
nificant fraction of GRBs (and, by association, of high- 
mass star formation) must occur within dusty regions 



not being probed by traditional optically-selected red- 
shift surveys. Based on our inferred distribution of Ay, 
we estimate that approximately ^50% of the rest-frame 
near-UV emission from Swift GRB afterglows is absorbed 
by dust. This value is quite similar to the fraction of ob- 
scured star formation inferred at high redshift based on 
far- infrared and millim eter studies (jCharv fc Elbadr2001l : 
iLe Floc'h et"all 120051 ) and may suggest that the poten- 
tial for GRBs t o serve as tracers of the hi gh-redshift star- 
formation rate (jBlain fc Natarai"anll2000( l is finally being 
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Fig. 9. — Rest-frame dust extinction Ay versus A'^^f as calcu- 
lated from absorption in the X-ray spectrum due to light metal ions 
(assuming Solar metallicity) for 24 of 29 bursts in the P60 sam- 
ple (five events, four dark and one non-dark, have been excluded 
due to the absence of meaningful constraints on either parameter). 
Bursts with known redshift are shown as solid points; bursts with 
unknown z are shown as open points at their most likely redshift 
(if only an upper limit is available, we plot the burst at a redshift 
of 2 ~ 2 or, in the case of GRB 080320, 2 ~ 5). A "track" hue then 
shows the evolution of the best-fit measurement or limit at differ- 
ent redshifts between z = 0.5 and the maximum host or afterglow 
redshift in Table [5] The majority of events have a ratio of Ay/Nn 
substantially lower than seen in Local Group galaxies, consistent 
with observations of other GRBs. (Milky Way, SMC, and LM C re- 
lations are plotted as lines using the values in lSchady et al.i |2007l. 
along with the average value for bursts in that paper and th e min- 
imum Ay/Njj in the pre-iSjui/t sample of IKann et al.|[2006l) . The 
high-Av events in our sample {Ay > 2), while not clearly in- 
consistent with the low Ay /Nu relation observed previously, may 
suggest a trend towards more 'normal' dust-to-gas ratios. 

realized. Nevertheless, there is still need for caution: in 
addition to the possibility that the extinction may be a 
unique property of the GRB site hinted at by the blue 
observed colors of the host galaxies in our observations, 
there is evidence that metallicity or other biases result 
in a GRB host population strongly f avoring subluminous 
galaxies in the local universe (e.g., Modiaz et al. n2008f) 
and possibly at highe r redshifts a s well [Fruchter et al.l 
[200a ILe Floc'h et all [2003; cf. iFvnbo et all 120081) — 
which could skew the GRB positional distribution signifi- 
cantly away from that of high-0 star formation in general. 

6. CONCLUSIONS 

Twelve years after the discove ry of a class of "d ark" 
GRBs lacking optical afterglows (|Groot et al.lll998[ ). we 
claim that the mystery surrounding the relative impor- 
tance of the varying hypothesis for their apparent optical 
faintness is largely resolved. Of 14 dark events (out of 29 
events in the full P60 sample): 

• Seven or more events (070521, 061222A, 060923A, 
060210, 080319C, 050416A, 080319A, plus likely 
050915A and perhaps 050713A) are significantly 
suppressed by dust extinction (at least 1 magni- 
tude in the observed i?-band and typically much 
more) in their z < 4 host galaxies. 



10.0 



0.1 







o Dcirk 
o Not dark 








^ : 







1 2 3 4 5 6 7 



z 

Fig. 10. — Constraints on rest-frame dust extinction (Ay) and 
redshift (z) inferred for all 29 bursts in the full P60 sample. As in 
Figure [9l ljursts with known redshift arc solid points; bursts with 
unknown z are shown as a redshift "track" showing the evolution 
of the best- fit Ay or limit with z; an open circle is plotted at a 
representative value. For a few bursts Ay is unconstrained; redshift 
limits are plotted as arrows at an arbitrary Ay with no circle. For 
clarity, the redshifts of two events have been adjusted slightly (less 
than 0.1) to prevent overlap of points. All bursts are constrained 
to 2 < 7 and all but one to 2 < 5 (for the exception, GRB 080320, 
extinction is not constrained above 2 ~ 6, as shown by the dotted 
line). However, many events show large extinction values, with a 
distribution skewed towar ds notic eably higher Ay than previous, 
nonuniform samples ("e.g.. iKann et al..,2007,y 

• Two events are probably suppressed due to Lyman 
absorption at redshift of z > 4.5: GRB 060510B at 
z = AMI and GRB 080320 (at unknown redshift, 
but z < 7). 

• Three events appear to be simply underluminous: 
not at high redshift, but because they were either 
intrinsically underenergetic (050607 and 060805A) 
or because little energy was coupled to the after- 
glow (050412, which may be a "naked" long GRB). 
Although no optical afterglow was detected for any 
of these e vents, they would not be classified as dark 
using the I Jakobsson et al.l (|2004a| ) criterion. 

We conclude that the dark burst phenomenon is pre- 
dominantly the result of extinction at moderate redshifts 
(1 < z < 5), with underluminous afterglows (other- 
wise normal events which are too faint for the sensitiv- 
ity of a small telescope) also contributing significantly 
in an amount depending on the detection threshold — 
consistent with, but more constraini ng than, the results 
of w e- Swift dark burst studies (e.g.. iDe Pa squale et al.l 
l2003f ). In particular, a large fraction of high-redshift 
GRBs is not needed, and in fact is ruled out — providing 
observational evidence limiting the ability of Population 
III stars to efficiently produce GRBs and in agreement 
with most recent models of the high-redshift star forma- 
tion rate. Furthermore, our methods suggest that even 
if the discovery of very high-z events continues to be 
extremely challenging (although the recent discovery of 
GRB 090423 at z = 8.2 has now demonstrated that such 
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events do exist and can be recognized), complementary 
host-galaxy searches can impose useful constraints on the 
high-redshift bursting rate free of selection biases, and 
we encourage continued host-galaxy follow-up of other 
medium-aperture robotic telescope samples and of dark 
bursts in general. (A much larger, though not uniformly 
sampled, broadband survey of Swift-era dark bursts is in 
progress.) 

The location and nature of this high-redshift dust re- 
mains unknown: although our wavelength coverage is 
limited, no galaxy in our sample shows evidence of sig- 
nificant dust extinction. In these cases, the line of sight 
to the afterglow must be passing through a much larger 
extinction column than the light from the observed young 
stars in the galaxy which dominate its rest-frame near- 
UV flux. The solution likely requires that the dust is 
nonuniformly distributed — either closely linked with the 
GRB site itself, or widespread but sufficiently patchy to 
conceal its effects. Although we cannot firmly resolve 
this question at this stage, it is clear that GRBs still 
have much to teach us about the structure of galaxies at 
high redshift and the importance of obscuration in the 
early universe. 

Fortunately, the answers to these lingering ques- 
tions may not be far off. Longer-wavelength observa- 
tions (near- and mid-infrared, sub-mm, radio) of these 
and other dark burst hosts would clarify the picture, 
piercing through the inferred dust screen or even de- 
tecting the reradiated emission from any postulated 
highly-extinguished population directly. Such stud- 
ies of the (liinited) pre- Swift daik bu rst host sample 
(|Barnard et al.ll2003l : iBerger et aLll2003f ) indicate a pop- 
ulation that differs little from GRB hosts in general — 
consistent with a patchy dust distribution in all GRB 
hosts, where the location of the GRB within its host 
(rather than the type of host) is the determining fac- 
tor in the observed darkness of a given burst. How- 
ever, the obscuration rates derived from these radio and 
sub-mm studies are surprisingly high (typical radio/mm- 
derived host SFRs are 20-50 times the optically in- 
ferred values) and the non-detection of most such sub- 
mm so urces in a recent surve y by the Spitzer Space Tele- 
scope (|Le Floc'h et al.ll200l ' may call this conclusion of 
ve ry high obscura t ion in to some doubt. The sample 
of iLe Floc'h eTall ((2006 1) includes three "dark" GRBs, 
one of which (GRB 970828) is indeed associated with a 
strongly obscured galaxy (one of only two identified in 
their sample of 16 objects). We suggest that more work 
in the long-wavelength regime is necessary to fully un- 
derstand the nature of GRB host galaxies, especially of 
the most highly-extinguished events. 

However, continued study in the optical band promises 
to be useful as well. High-resolution space-based imag- 
ing could constrain the morphologies of dark GRB host 
galaxies, including any possible difference between the 
burst site and the rest of the galaxy. (For example, 
studies of the host of pre-Swift dark GRB 970828 seem 
to indica te a dust lane running through the afterglow 
position [Diorgovski et al.l 1200 1| .') The most luminous 
bursts are capable of shining through even very thick 
dust columns, allowing for detailed study of the mate- 
rial along their lines of sight. Recently, spectroscopy 
and infrared photometry of the extremely bright GRB 
080607 at z = 3.036 revealed a strongly extinguished 



(Ay = 3.2 ± 0.5) event, showing a clear 2175A bump 
and an abundance of molecular and ionic lines associ- 
ated w ith a nearby molecular cloud with Solar- like metal- 
licity (jPro chaska et al. 2009). Similarly, S'wi/i( bursts 
GRB 050401 (Wa tson et all 12006 ) and GRB 070802 
(Kriihler et al.. . 20081 ) were also "dark" events that were 
nevertheless sufficiently optically luminous to enable 
multiband photometry and optical spectroscopy, con- 
firming the link between optical suppression (darkness), 
reddening, and dust absorption. Such events, while rare, 
illustrate the need for continued observational effort on 
as many fronts as possible (including both spectroscopy 
and photometry, of both afterglow and host galaxies, 
and at all available wavelength regimes) to make further 
progress on the environments of gamma-ray bursts and 
their connection to star formation in the early universe. 
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